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PLASTIC DEFORMATION OF WIDE- FLANGE BEAM-COLUMNS 


Robert L. Ketter,* Edmund L. Kaminsky** and Lynn S. Beedle*** 


I. Synopsis 


Presented in this paper are the results of an analytical and experimental 
study of the inelastic deformations of wide-flange steel beam-columns. 
Emphasis has been directed toward the influence of axial thrust on the 
moment-curvature (M-@) relationship when stresses exceed the elastic limit. 
The analysis is extended to include the influence of residual stresses, and 
thus the moment-curvature relationship is developed as a function of axial 
thrust and residual stress. Buckling loads are then determined for eccen- 
trically loaded steel columns of WF shape, the necessity for an assumed 
‘shape factor’ being eliminated. These derived moment-curvature relation- 
ships are also applied to determine column end rotations and mid-height 
deflections. 

Two sizes of WF shapes are considered, (BWF31, 4WF13) and in the ex- 
perimental studies, tests were carried out on as-delivered WF structural 
sections. 


II. Introduction 


When extending the methods of analysis used in engineering mechanics to 
include the case of elastic-plastic bending, the relation between applied 
moments, axial thrust and resulting curvature for the case where stresses 
exceed the elastic limit is basic in importance, since it is the integrated 
effect of this curvature that determines deflections and rotations. The value 
of moment-deflection data is primarily in that it allows the determination of 
true maximum loads and the deflection at those loads. Moment-rotation data 
is also valuable since it enables one to evaluate performance relating to 
energy absorption. Further, it indicates the amount of rotation that a member 
may sustain at near-maximum moments. 

For flexure in the elastic range the curvature, ¢ , is given by 


=M/EI (1) 
where 
¢ = curvature, 
M= bending moment at the section in question, 
E = Young’s modulus of elasticity, and 
I = moment of inertia of cross-section about axis 
perpendicular to plane of applied moments. 


Equation (1) can be extended (neglecting certain higher order terms) in the 
form of the differential equation 
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By direct integration, Equation (2) defines all the various load-deformation 
relations. 

On the other hand, above the elastic limit such ease of solution does not 
exist, since stress is no longer a linear function of strain and, therefore, od 
does not vary linearily with M. The problem of determining M-@ relations 
in the plastic range has been solved for the condition of pure bending. 10 
Since for that problem the neutral axis of strain coincides with the geo- 
metrical center of the section, relatively simple expressions were developed 
relating M and¢ . If, however, axial load is applied in addition to the im- 
posed moment, such a strain condition no longer exists. 


1. Influence of Axial Thrust 


The determination of M-¢ curves including the effect of axial load has 
been developed by Timoshenko for the case of rectangular sections.* For 
the case of WF shapes, discontinuities of the cross-section make direct com- 
putation troublesome, and in the next section of this report a ‘point-by-point’ 
method of computing M-¢ curves for rolled shapes is presented. It is similar 
to the method used earlier for the pure bending case,!0 but includes also the 
influence of axial thrust. 

An important assumption made in this derivation is that the Bernoulli- 
Navier hypothesis (bending strain is proportional to the distance from the 
neutral axis) can be extended to include the case of plastic deformations. ** 
Using this hypothesis, together with assumptions with regard to properties 
of the material in question, this report outlines and illustrates a method 
whereby a relation between axial load, bending moment and curvature can be 
obtained even though stresses exceed the elastic limit. 

The curvature, @ , is a function of the strain developed at a section; the 
moment, M , and the axial thrust, P , are functions of the resultant stress. 
Fig. 1. Then, 


P- odA (3) 


M = Be (4) 
A 


@= (5) 
d 
For the case of elastic strains illustrated in Fig. 1 (that is, when €, re 


and € 9 <4 y the distribution of strain and distribution of stress are 


* Pages 51-54, Ref. (12) 

** Actually for steel members this is an idealization. As has been described 
in previous work (16), it is only required that in the elastic range, the 
strain be proportional to the distance from the neutral axis, and in the 
plastic range, that the stress be equal to the lower yield-point stress. 
Thus the curvature at any section is a function of the part which remains 
elastic. This is further discussed in Section III of this report. 


If the strain-hardening range is considered, it is once again necessary to 
make an assumption with regard to strain distribution in the inelastic 
range. But strain-hardening has not been considered in this paper. 
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Similar, 9 1 being equal to E €4, etc. When yielding occurs, however, this 


linear relation no longer holds; instead, there results a stress distribution 
composed of segments of various straight lines or curves depending on the 
stress-strain properties of the material in question. 

In the analysis the idealized stress-strain diagram of Fig. 2 is assumed. 
This curve closely approximates the experimental stress-strain diagram of 
mild structural steel in the elastic and plastic ranges but strain hardening 
is neglected. Figure 3 is typical of the resultant stress and strain distri- 
bution curves when strains exceed the elastic limit. The stress is equal to 
the yield-point value when € > € y’ As would be expected, the moment 


corresponding to a given value of ¢ would be less than that predicted using 
the straight-line relation of the elastic solution. However, by assuming a 
given stress pattern, it is possible to compute directly the corresponding 
value of @ ; likewise, P and M may be determined. This may be done for 
various stress patterns and a curve plotted relating the three variables, P, 
M, and @ .* 


2. Residual Stresses 


The behavior of structural members may be modified considerably by the 
presence of residual (‘locked-in’) stresses. A short review of their nature 
and influence is now presented. 

Residual stresses are of two basic types: those due to differential cooling 
of the member during and after the rolling process, and those due to cold 
bending of the specimen while being straightened or during fabrication. 


Cooling Residual Stress 


When the ingot is originally placed in the rolling mill, the temperature is 
nearly constant throughout the material. However, as it traveis through the 
various rolls (each reducing the thickness of the different parts of the 
section), and later on the cooling bed, a differential in temperature develops. 
As a result of this differential and the subsequent plastic deformation of dif- 
ferent parts of the cross-section at high temperature, the parts that cool the 
fastest are usually found to be in residual compression while the slower cool- 
ing portions indicate tensile residuals. 

For a flat plate, the resulting stress distribution pattern after thermal 
equilibrium might be of the form shown in Fig. 4. Extending this reasoning 
to a rolled WF shape, the cooling residual stress pattern might be similar to 
that shown in Fig. 5. The mechanism, by which these residual stresses are 
formed has been described in several references.“ The influence of 
cooling residual stresses on axially-loaded specimens has been described in 
a recent paper.6 

The existence of such initial internal cooling stresses in ‘as-rooled’ stee! 
members is demonstrated in several common shop practices. One is that of 
cutting a WF section along the center of the web to produce two T sections. 
As the cut progresses from one end of the beam toward the other, the indi- 
vidual T’s deform--curving away from or toward each other, depending upon 
the distribution of residual stresses. 


* For the rectangular section it is not necessary to resort to graphical repre- 
sentation of the inter-relation between Equations (3), (4) and (5) since ana- 
lytic expressions connecting these three variables are readily available. 
(See, for example, equations (42) and (47), page 37 and 39 of reference 4 or 
reference 13, page __i+). 
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As a basis for making an assumption with regard to the magnitude and 
distribution of residual stress, Fig. 6 shows a typical plot of measured resi- 
duals in the flanges and web for one of the sections considered in this investi- 
gation. The experimental points could be approximated either by parabolas 
or by straight lines. The latter has been chosen to facilitate computations 
and current work shows that this is a reasonable approximation for the shape 
studied. 6 

As re-plotted in Fig. 7, across each flange it is assumed that residual 
stress decreases linearly from a maximum compressive stress, g ae? at 


the flange edge to a tensile residual stress, at the web junction. 


st’ 

Furthermore, it is assumed that this tensile stress is constant across the 

web. It is required that og oe and g ot be chosen in such proportion that 


static equilibrium of the cross-section is maintained. This condition is as- 


sured if, 


rt bt R, 


Tre bt*wid-2t) (5.1) 


b = Flange width 
t = Flange thickness 
w = Web thickness 

= Section depth 


Residual Stresses Due to Cold Bending 


After the member passes through the last rolling mill and is cut into 
lengths suitable for handling, it is placed on the cooling bed. During cooling, 
the member tends to deform due to uneven support and other causes and, 
when cold, is no longer straight. Therefore, to obtain a straight finished 
product, each member must be cold bent. The conventional straightening 
method is to place the member in a ‘gag’ and plastically deform the curved 
part in the opposite direction to the original curvature. This ‘gag’ is dia- 
grammatically shown in Fig. 8. 

A member enters the ‘gag’ bent as shown by the solid line in Fig. 8. By 
applying a load, it is forced to the position indicated by the dashed lines. If 
on release of this load, the member is to have the straight form shown by the 
dash-dot position, a certain amount of yielding must take place to enable the 
‘permanent-set’ to compensate for the original curvature. Since yielding has 
taken place, the strains in certain of the fibers have exceeded the elastic limit 
of the material and will, therefore, not return to their original positions on 
removal of the load. The balancing of these strains across the section pro- 
duces cold-bending residual strains. Progress Report No, 515 contains a 
discussion of the magnitude and distribution of residual stresses formed as a 
result of cold-bending. Fig. 9 is typical of a ‘cold bend’ residual stress 
pattern for rectangular members (WF sections deformed about their weak 
axis). 

Since plastic deformation is localized at points of load-application in the 
‘gag’, it would appear that cold-bending residual stresses are also local in 
nature. At any rate, it is considered that these initial stresses do not have 
as pronounced an effect on column behavior as do cooling stresses which are 
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distributed rather uniformly along the length of the member. Therefore, de- 
tailed consideration of the influence of cold-bending residual stresses will 

not be made, emphasis being directed toward the cooling type of initial ‘locked- 
in’ stresses. 


3. Stress Distribution Beyond the Elastic Limit 

While the determination of stress-distribution and resulting moments and 
loads from an assumed yield condition is straight-forward for members with- 
out residual stress, the superposition of load-carrying strains and residual 
strains requires special attention for members that do contain residuals. 

Consider the rectangular section of Fig. 10a which contains a certain resid- 
ual strain (and stress) pattern. To this pattern of initial strain superimpose 
the strain pattern of an externally applied bending moment and thrust (defined 
by the applied strain pattern of Fig. 10b). If the total strains resulting from 
this superposition (as shown at the top of Fig. 10c) are within the elastic limit, 
the applied stress and applied strain patterns are similar across the entire 
section. Since @ is a function of strain whereas M is a function of stress, 
it is clear that Equation 1 is applicable to this bending problem. When the 
combined strains exceed the elastic limit (Fig. 11lc at top) @ is no longer a 
linear function of M. Similarity between stress-and strain-distribution exists 
only for those regions in which the combined strains are less than the yield 
point (shown as a in Fig. 11b). The resulting stress pattern (Fig. 11c) is 


determined directly from the strain pattern above. The applied stress pattern 
(Fig. 11b) is obtained by subtracting the initial residual stress pattern from 
the ‘resulting stress pattern’. (This description is applicable to a WF column 
bent about the weak axis. . . .so long as the web is neglected.) 

Extending this reasoning to the case of flexure of a WF shape bent about its 
strong axis, the result of addition of residual and load-carrying strains is 
illustrated in Fig. 12. On the basis that plane sections before bending remain 
plane after bending, the applied stress patterns shown in Fig. 12 may be de- 
termined by the subtraction process just described. Under the applied moment 
and axial load shown diagrammatically, the stress-distribution for Case 1 is 
at the instant of first yielding of the compression flange. In Case 2, this 
flange is partially plastic and the tension flange has just reached the residual 
yield stress. When the compression flange is completely plastic (Case 3), 
the tension flange is still partially plastic. In Case 4, the member is com- 
pletely plastic.* 

Similarly, for the cold-bending type of residual stress (cold bending as- 
sumed to have been about the weak axis of the section), ‘stress addition’ is 
shown in Fig. 13. Again, linear superposition of strains is assumed. 

It is the purpose of this report to present a method for computing the 
moment-curvature relationship for WF columns, bent about either axis, and 
in the presence of axial thrust. It will also be shown how the presence of 
residual stress modifies the M- @ relationship. Correlation with tests will 
be shown and the results will be used to predict the strength of as-delivered 
columns loaded with equal eccentricities. In the theory, the possibility of 
local and lateral-torsional instability is not considered. 


* P-M relationship is as described in Progress Report No.6(9) 
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Ill. M- Relationship 


1. Members without Residual Stress 
When a member is loaded within the elastic range by combined bending 
and thrust, stress and strain diagrams similar to those shown in Fig. 14 
result. Since the strains are small in comparision to the depth of the section, 
1- €2 01. %3 


For members loaded beyond the elastic limit (Fig. 15) 
1. *3 €y - 


oy 


It is thus possible to compute ¢ directly from a stress diagram by consider- 
ing only that part of the section which remains elastic. 

From a stress diagram, it is also possible to determine both the moment 
acting on the section and the applied thrust.!2 By using Equations (3), (4) 
and (5b) it is then possible to determine M, P and ¢ for any given stress dis- 
tribution. 

As loads are increased under a given loading condition, it is likely that 
successive stress-distributions will be similar to those shown in Fig. 16. 
Yielding will first occur at the extreme fiber on the compression side of the 
specimen. As the loading is increased, yielding will penetrate through the 
flange on the compression side. Eventually, yielding occurs on the tension 
side--Sketch (c)--and the full plastic condition is typified by Sketch (d). 

For ease of computation certain flexure conditions are assumed, and a 
typical set of these is illustrated in Fig. 17. From the equilibrium conditions, 
values of P and M may be determined for each assumed value of g@. This 
makes it possible to plot curves of P vs ¢ and M vs @ for any given value of 
yield penetration on the compression side. One such set of curves is shown 
in Fig. 18 for an 8WF31 section where yielding was assumed to have pene- 
trated to a depth of 1/4d from one side. (It is Case (c) of Fig. 17.) These 
curves have been made non-dimensional by plotting P/P, vs $/dy and M/My 
vs@/%,. For all points plotted, the corresponding stress distribution dia- 
grams have been shown. 

Curves similar to Fig. 18 have been computed for each of these assumed 
values of yield penetration for the 8WF31 section and are presented in Fig. 19. 
Since the straight-line portion of each of these curves is that corresponding 
to an elastic stress condition on the tension (convex) side of the section, only 
the limiting cases where 9 g = yield point stress in compression and tension 
need be investigated to define that part of each curve. The limit of this range 
has been shown as a dash-dot line. 

The method of using these ‘auxiliary curves’ (Fig. 19) is to select a value 
of P/Py for which a M-@curve is desired. For example P = 0.2P,. With 
this vales of P/Py draw a horizontal line cutting the various auxilfary p-¢ 
curves. At these points of intersection, project verticals until they intersect 
the corresponding M-@ auxiliary diagrams. Connecting these points gives 
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the desired M- @ curve including the effect of an axial thrust of 0.2Py. (This 
process is indicated in Fig. 19 by dashed lines.) Fig. 20 is the desired curve 
with the auxiliary construction lines removed. 

The curved portions of Fig. 19 have been plotted to a larger scale in Fig. 
21 to facilitate greater accuracy in this graphical miethod. Here M-¢ curves 
for several different values of P/Py have been shown (dashed lines). 

The influence of varying axial load for a range of P Py is illustrated in 
Fig. 22. M-@ curves have been plotted here for various values of P/P 
ranging from 0 to 0.8. As indicated by the dashed line, only when P is rela- 
tively small will the tension flange be plastic. For the range in which it re- 
mains elastic, only the straight-line portion of each auxiliary curve need be 
considered to obtain the desired M- @ curve. 

Table 1 contains equations necessary for plotting critical points on the 
‘Auxiliary curves’ (Fig. 19 or 21). It will be recognized that the computation 
is time-consuming and tabular solution appropriate. 

No tests using annealed material were conducted as part of this investi- 
gation and consequently comparison between the theory and test is not made 
here. However control tests on axially-loaded columns were carried out on 
the 8WF31 material as part of another study®. 


2. Members Containing Residual Stress 

By a process analagous to that just described, expressions for P, M and @ 
may be developed for assumed ‘yield conditions’ for sections containing resid- 
ual stresses. (See Fig. 12 for flexure about the strong axis, Fig. 11 for weak 
axis bending). The equations for flexure about the strong axis are shown in 
Table 2 and those for flexure about the weak axis are given in Table 3*. In 
these equations, 


= Fre /%y 
= Trt /%y 
* 
/a 


tt 
The equations summarized in Table 2 reduce to those shown in Table 1 for 
annealed material when the residual stress is zero (R, = R = 0). 


Auxiliary P-¢@ and M- ¢curves are shown in Fig. 23 for one ‘alpha’ con- 
dition (a@ = 0.4 indicates that yielding has penetrated to four-tenths the depth 
of the member). The assumed residual stress function R,, = 0.3) is con- 
sistent with actual residual stress measurements. The Figure also shows 
the auxiliary curves for material free from residual stress. These curves 
have also been made non-dimensional by plotting P/Py - vs- $/Py and M/My 
- vs- Y Dy- The portions which are cross-hatched correspond to regions 
covered by the equations of Tables 1 & 2. Where shading is missing, the 
member has partially yielded in the tension flange and the equations for this 
short region are not presented. 

This method of assuming a penetration of yielding on the compression side 
of the section and then determining curves of P -vs- g@ and M -vs- ? can be 
carried out for any value of yield penetration. Six different values were 


* In the derivation of these weak axis equations the web was neglected and the 
section was considered as two equivalent rectangles. In reality then these 
expressions of M, P and ¢ apply strictly to the rectangular section. 
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assumed to obtain the curves presented in this report. For each, a set of 
curves similar to those of Fig. 23 was plotted preparatory to plotting the M-? 
curves. The latter are obtained for a given axial load by drawing a horizontal 
line cutting the various auxiliary P-% curves at the assumed axial load value 
(P/Py = 0.4 has been chosen for illustration). At these points of intersection, 
project verticals until they intersect the corresponding M-® auxiliary dia- 
grams (circled in Fig. 23). Plotting the values obtained from each of the 
separate ‘auxiliary curves’ and connecting the corresponding points gives the 
desired curves. These are plotted in Fig. 24, the points obtained from Fig. 23 
being circled for illustration. 

Typical M- @curves for flexure about strong and weak axes for the 8WF31 
section are shown in Figs. 25 and 26, respectively. In each of these figures, 
solutions neglecting the influence of residual stress are shown by solid lines 
whereas those including this variable are indicated by dashed lines. As ex- 
pected, the influence of residual stress on M- ¢ curves is more pronounced 
for flexure about the weak axis. 


3. Measurements 

In Figures 27-31 are plotted the M- % relationships obtained in the course 
of testing sixteen different columns and Table 4 summarizes the pertinent 
loading information for each of the tests reported. It also indicates for each 
specimen the observed mode of failure. All specimens were tested in the 
as-delivered condition and two sizes were used (8WF31 and 4WF13) in a range 
of L/r from 21 to 111. Details of test apparatus and procedure are described 
in Ref. 3 and Ref. 9 summarized the behavior of the specimens tested. 

The curvature, @, at any section is a measure of the variation in strain 
differences across the section. SR-4 strain gages were placed symmetrically 
about the centerline of the section at a known distance from one another, this 
distance being such that 6-measurements could be made relatively far into 
the inelastic range. They were placed as close as practicable to the section 
at which maximum moment was expected. 

Moment values at ‘¢-sections’ were determined as the sum of applied end 
moments and the product of applied axial thrust and measured deflections. 

In Figs. 27 through 31 both experimental and analytical M- @curves are 
presented.* With each set of curves is shown the loading condition with the 
column in its deflected configuration. The shaded arrows denote the location 
of the section at which @ was measured. The other auxiliary sketch is an 
interaction curve? which shows the maximum values of axial thrust and end 
bending moment at collapse as well as the path followed in obtaining that 
loading. 

The theoretical M-@ solutions neglecting residual stress are shown by 
dashed lines. For an assumed flange tip residual stress equal to three-tenths 
of the yield stress (i.e., R, = 0.3), the theoretical curve is shown by a dash- 
dot line. The experimental curves are shown solid. Initial yield as determined 
by flaking of the mill scale coated with whitewash is shown for each test by the 
open arrows. 

In the figures the symbol ‘LA’ denotes lateral buckling as observed with 
the aid of dial gages and SR-4 gages; the symbol ‘L’ indicates oberved local 
buckling. Point ‘LA’ was well defined by the measurements made. (The mode 


* The material properties for the two sections investigated are tabulated in 
Table 5. 
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of failure shown in Table 4 is an effort to show how the column finally col- 
lapsed). In each case, lateral buckling occurred at a moment higher than the 
first observed inelastic action, but the onset of lateral buckling is in most 
cases not the maximum observed moment. Rather, it is a point at which de- 
formations commence to increase at a more rapid rate. For most of the 
specimens, the test results are in fair agreement with the theory up to the 
point at which lateral-torsional buckling was observed. 

Referring to Fig. 28, the location of the open arrows suggests that the 
residual stress level was about the same in each of these specimens (Re = 
0.26). This is interesting because two different rollings from the same ingot 
were involved. Further, since lateral buckling (‘LA’) was observed for each 
of these condition ‘c’ tests at about the same moment and since the slender- 
ness ratios are different (L/r varies from 27 to 55), inelastic lateral 
buckling in this range appears to be a function of maximum moment and is 
not much affected by the length. 

The behavior of T-13 (Fig. 29) is different from the rest of the eight-in. 
specimens, the curve closely approaching that for a member without residual 
stress. A possible explanation is that the member might have been cold- 
bent in the region at which @ was measured thus wiping out the cooling resid- 
ual stresses. 

In Fig. 31, showing the results of tests on four-in. members, the open 
arrows indicate that this material had a lower level of residual stress than 
did the eight-in. material. (The average value for all 4-in. specimens tested 
was R, = 0.1 whereas that for the 8-in. specimens was R, = 0.24). 

Specimens T-3 (Fig. 27) and T-21 (Fig. 30) were identical except for the 
size of member. The measured M- ¢curves are also similar, deviating 
markedly from the theoretical curve. This seems to be characteristic be- 
havior when a relatively greater axial load is applied to the member. Also 
owing to the fact that the column is effectively restrained by the condition 
‘b’ loading, lateral buckling is not observed until a considerable amount of 
plastic deformation has occurred. Incidentally, the reduction in strength of 
the theoretical curves to T-21 below those of T-17 and T-20 indicates the 
influence of axial thrust summarized earlier in Fig. 25. 

Good correlation between theory and test is observed in T-7 and T-9 (Fig. 
29) for 4-in. material even though two different rollings from the same ingot 
are involved. 

In Fig. 32 are plotted M- curves for several different tests. It is seen 
that the experimental results are consistent. 

Examining the individual curves for tests 12, 14, 16, and 19 (Figs. 27, 28) 
the M- ¢@ behavior is relatively consistent for these 8WF31 specimens below 
the point at which twisting of the specimen was noted. It would not be expected 
that the theory would apply beyond the point at which buckling failure occurred. 
However, even prior to lateral-torsional buckling the curves did not corre- 
spond exactly to theory. Several reasons for this discrepancy are: 

(a) The assumed residual stress pattern ‘straight-line variation) is an 
idealization. A variation in magnitude of residual stress also will affect the 
M- ¢curve. 

(b) It was assumed that yielding progressed uniformly to known, given 
depths (Ycec and yy); however, steel does not behave in this manner, certain 
parts yielding more and others yielding less than ew This problem has 
been extensively discussed in Progress Report No. 5 5, and Fig. 33 selected 
from that report suggests that curvature for a given moment would be greater 
than predicted by the theory described earlier in this paper. 
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IV. Column Strength 


The results obtained in the previous section will be applied to a consider- 
ation of the eccentrically loaded member shown in Fig. 34. This single curva- 
ture loading condition is essentially the same as loading condition ‘c’ (Table 
4). The member is assumed to be initially straight, of uniform cross-section 
and free from accidental eccentricities. For WF shapes it is of interest to 
determine the following: 

(a) the ultimate strength of the eccentric column 

(b) the load-deflection relationship for columns with equal eccentricities 

(c) the end-rotation of columns in the plastic range. 

Knowledge of the M- g relationship makes possible a solution to these prob- 
lems--including also the influence of residual stresses. 

The strength of the eccentrically-loaded column has been considered by 
Bleich, who reviewed the work done up to that time.4 As he pointed out, a 
determination of true column strength depends on a knowledge of the load- 
deflection relationship (Problem ‘b’, above). End rotations (problem ‘c’) are 
possibly of academic interest for columns loaded with equl eccentricities. 
Such members collapse quickly after reaching the maximum load. They have 
little ‘rotation capacity’, i.e., the ability to deform through considerable angle 
changes at near-maximum loads. On the other hand, columns bent in double 
curvature may show considerable ‘rotation capacity’, and indeed this is re- 
quired if structural members are to meet the assumption implicit in plastic 
analysis and design.* 


1. Load-Deflection Relationship 
The difference between the elastic and the plastic load-deflection re- 
lationship is shown in Fig. 36. The equation of the elastic curve is 


A __i1-cosu 
e  cosu 


where 
u = kL/2, e = eccentricity of load, A = centerline 


deflection, and k = YP/EI. 

Even with large eccentricities, the maximum elastic strength approaches 
the Euler buckling load. However, as the load increases, the stress in the 
outermost fibers at the center section of the column increases at a faster 
rate due to the bending moment which is increasing as the product of P and 
the quantity (A +e). See Fig. 34. Therefore, columns of the proportions 
used in engineering structures start to yield, say at some point ‘a’ in Fig. 36. 
As the compressive force, P, is increased beyond this initial yield value, 
plastic flow spreads from the edges towards the center of the section and 
along the bar thereby reducing its flexural rigidity. Consequently, the de- 
flection, y, (and hence the moment of the external forces) increases at a 
' higher rate than for the elastic case. Finally a point is reached at which 
further increase in load is impossible and collapse occurs. This maximum 
load would correspond to point ‘b’, the portion ‘b-c’ representing the inelastic, 
unstable load-delfection relationship. 

The load deflection curve is determined by integration of d-values. How- 
ever, ¢ is dependent upon the bending moment which in turn is dependent 


* This is discussed in several papers, most recently in Ref. 8. Column end 
rotations have been examined in a recent unpublished report, ‘Moment- 
Rotation Characteristics of Beam-Columns’, by R.L. Ketter and L.S. Beedle, 
Progress Report R, Fritz Laboratory Report 205A.11, November, 1952. 
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itself upon the deflection. A numerical integration procedure!} is suitable 
in problems of this type. A deflected shape is assumed and moments are 
computed. ¢ -values are obtained (Section III) and these are integrated to 
obtain a new deflected shape. If the assumed curve is correct, the deflection 
values obtained will be identical with the former. This method was applied 
to the members studied in the program and Fig. 35 shows typical resulting 
curves. The solid curve is the experimental result and the theoretical curves 
using the numerical method are shown by dot-dash lines for zero residual 
stress and for residuals estimated at three-tenths of yield-point stress. 
Agreement is reasonable up to the point at which the member starts to twist. 
Even though the exact deflected configuration is that which results in the 
lowest potential energy of the system, the magnitude of the critical load is 
relatively insensitive to changes in the assumed deflection curve providing 
the deflection patterns conform to the geometric boundary conditions of the 
member. Bleich? describes the use of an assumed sine curve by Jezek and 
an assumed partial cosine curve by Osgood and Westergaard. These two 
assumptions have also been used in the following paragraph toward the de- 
termination of column buckling loads in the inelastic range. 
Assuming that the column bends in a sine curve, 
x 
y= A sin = 
The section at mid-height is subjected to the greatest deformation, and ¢ 
at this section is d2y 
(7) 


2 
dx L 


The moment at this same section for the eccentrically loaded column is given 
by 
M=P (A+e) 


and for the condition- ‘c’ -loaded column by 
(9) 


From Eq’s. (7), (8) or (9) and the M-¢ plot for the section in question, load- 
deformation curves may be obtained. A typical resulting curve is shown in 
Fig. 35 by dotted lines for the two material conditions (annealed and as- 
delivered). The predicted maximum end moment agrees with the numerical 
method within one percent. 

Referring to Fig. 34, the deflection expressions for a partial cosine curve 
is given by 


(8) 


4 
y=(A4 +e) (10) 


Then, by using the first two terms of a cosine series, 


(11) 
L 


In this case the resulting end-moment vs deflection curve coincides almost 
exactly with the numerical integration curve of Fig. 35 and has therefore not 
been shown. 

Reviewing the detailed steps which were followed in obtaining these ap- 
proximate load-deflection curves; 

(a) Assume a center line deflection, A . 

(b) Using either equation (7) or (11), compute gd. 
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7 (c) Enter the appropriate M-¢ curves with the value of o/b, and read off 
the corresponding M/M, ratio. 
(d) Compute M,, or P from either equation (8) or (9). 
2. Strength of Eccentric Columns 
: The ordinary column curve presents strength as a function of slenderness 
ratio, L/r. By using the sine curve assumption and the M-¢ curves of Fig. 
25 and 26, column curves are obtained as presented in Fig. 37 (strong axis 
Shox bending) and Fig. 38 (weak axis bending). Data is shown on a non-dimensional 


basis, the maximum load being divided by Py determined from the product 
of yield-point stress and cross-sectional area. The curves show the influ- 
ence of end eccentricity and of residual stresses. Three curves are shown 
for each eccentricity ratio (ec/r“): the maximum strength of annealed 
members, the maximum strength of as-delivered columns containing residual 
stresses and the initial yield (secant) solution. The shaded portion at each 
eccentricity shows the influence of residual stress on collapse strength. 
Curves of allowable load according to the A.A.S.H.O. Specification! are also 
shown. 
The curves show that the secant solution is on the unsafe side for inter- 

: mediate length columns with small eccentricities and that it is too conserva- 

ah tive for short members, the latter being particularly true for large eccen- 
. tricities. For small eccentricities, columns will collapse at loads consider- 

5 ably less than predicted by the initial yield condition. At zero eccentricity 

. a the considerable reduction due to residual stress has been confirmed by 
tests and is reported upon separately.6 

Fig. 39 is a cross-curve of Fig. 37 (strong axis) at L/r=55. In addition to 

the various theoretical solutions the results of three tests are shown. The 
collapse strength for each of these test members is less than the theoretical 
solution that includes the influence of residual stress. This is due, in part 
at least, to flexural-torsional buckling of columns bent about the major axis, 
ie / the mode of failure observed for each of these tests. 

The ‘collapse’ curves of Fig. 37, 38 and 39 are based on the assumption 
that the section deforms only in the plane cf the applied moments and that 
during this deformation it maintains its original cross-sectional form. While 

i this is satisfactory for weak axis bending, lateral-torsional buckling occurs 
i j when flexure is forced about the strong axis of the section. A solution to this 
problem is to be sought in the program currently underway at Lehigh. 
: An opportunity to compare the theory with the results of tests on annealed 
a si columns is afforded by experiments carried out recently at Cornell!4. In that 
‘g investigation columns were allowed to bend about the weak axis. Fig. 41 shows 
: 4 the comparison between the theoretical solution outlined above and the re- 


ta Ps sults of experiment (circles). The theoretical curve is obtained by making a 
a cross-curve of Fig. 28 at an L/r of 50, the I-shape being approximated by 
Pr two equivalent rectangles. It is seen that the agreement is good. As-delivered 


specimens were also tested, but since the level of residual stress is not 
known, the corresponding test points and the theoretical curve have not been 
plotted, although a rough approximation to a possible residual stress pattern 
yields a fair correlation between test and theory. 


3. Moment-Rotation Relationship 


_ Typical end-moment vs end-rotation curves are shown in Fig. 40. The 
om R theoretical curves obtained by the method of numerical integration are shown 
* by dashed lines, the test result shown by the solid line. For the theory that 
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takes into account residual stresses, and in the region prior to twisting of 
the cross-section, agreement with test is excellent. 


V. Discussion and Summary 


1. M-@ Relationship with Axial Thrust Present 

A method has been presented for determining the influence of axial load 
on the moment-curvature relationship in the plastic range. Simultaneous 
thrust-, moment-, and curvature-values are determined (Table 1) for dif- 
ferent yield conditions, and from these are constructed ‘auxiliary curves’ 
(Fig. 19). Once these auxiliary curves are plotted for a particular wide- 
flange shape, the M- ¢ relation may be obtained conveniently for any desired 
value of axial load (Fig. 22). 

A significant reduction in moment-carrying capacity at a given cross- 
section is to be expected in the presence of axial thrust. (Fig. 22). 

Previous solutions were limited to rectangular shapes. This graphical 
method extends the solution to Wide-flange and other symmetrical shapes. 
Since all rolled shapes contain some residuals, the solution has practical 
application in that it may be extended to members containing such ‘locked- 
in’ stresses. 


2. Influence of Residual Stress 

The above method may be applied to members containing residual stresses 
after a distribution pattern for such stresses has been assumed (Fig. 6, 7). 
For flexure about the strong and weak axis, equations are developed (Tables 
2 and 3) for the construction of auxiliary curves (Fig. 23) from which M- ? 
curves can be obtained by graphical construction (Figs. 25, 26). 

For a given axial thrust, residual stresses increase the deformation in 
such a way that it appears as though the moment strength were reduced. If 
straining can be continued sufficiently without column buckling occuring, the 
effect of residual stress is wiped out and the strength reaches the predicted 
value (Fig. 25). Tests show that column collapse may occur in the region in 
which residual stress is of influence (Fig. 28), and this is discussed further 
in the next section. 

The influence of residual stress is more pronounced for members bent 
about the weak axis than for those bent about the strong axis (Fig. 26 com- 
pared with Fig. 25). 

Test results in general confirm the theory up to the point at which lateral 
deformation commences for members bent about the strong axis. ‘Weak- 
axis’ experimental data will be available later in the program which is 
currently underway. 


3. Column Behavior 

Application of the M.-¢ data derived in this report to compute load- 
deflection curves for eccentrically-loaded columns (by numerical integration) 
gives reasonable correlation with test results up to the point of twisting 
(Fig. 35). A nearly identical solution is obtained by assuming a partial cosine 
curve. 

Excellent correlation between theory and test (up to the point of twisting) 
is obtained for the relationship of end bending moment to end rotation of the 
column (Fig. 40). 

The ultimate strength of eccentrically loaded columns may be determined 
by the maximum point of a load-deflection curve. The results obtained using 
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such a procedure are summarized in the column curves of Figs. 37 and 38 
and show the maximum strength of annealed members, of members containing 
residuals, and the yield stress (secant) solution. The reduction in collapse 
strength due to residual stress (shaded) at L/r = 85 ranges from 22% for 
ec/r2 = 0 to about 4% for each ec/r2 = 1.0 for the strong axis. Similarly 
columns which collapse due to bending about the weak axis suffer theoretical 
strength reductions of 26% at zero eccentricity to 6% at ec/r2 = 1.0. 

The initial yield solution (secant) is generally on the unsafe side for inter- 
mediate length members with small eccentricities, and is too conservative 
for short members. 

Fair correlation between theory and test is established in Fig. 39 fora 
few columns bent about the strong axis and which eventually collapsed by 
lateral-torsional buckling. 

Final failure for all members tested in the entire program was due either 
to lateral-torsional buckling or to local buckling (Table 4). However all of 
the theoretical curves derived and presented herein assume that failure is 
due to simple column bending. A consideration of test results (see Fig. 40, 
for example) suggests that the influence of residual stress might be more 
than the initiation of yielding at low loads with the subsequent loss in buckling 
strength shown in the theoretical curves. This early yielding due to residual 
stresses may aggravate the tendency toward lateral and/or local buckling; 
further studies are being made on the problem of inelastic instability. 

A considerable amount of computation and plotting is involved in obtaining 
these results; simplified expressions are required which will sufficiently de- 
scribe the behavior of these columns. Work has been started to apply the 
results to columns under different loading conditions and with end restraints. 
The column strengths here reported represent a ‘worst case’. Higher 
strengths will be expected for other loading conditions and end restraints. | 

Whatever column analysis is made, the basic relationship is that of 
moment to curvature. The results presented here are dependent upon an 
assumed magnitude and distribution of residual stress. The assumption seems 
reasonable, in the light of measurements that have been made up to this time. 
Eventually, when more data is collected, current design rules may be re- 
examined toward providing more economical structures consistent with 
adequate safety. 
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VII. Nomenclature 


Area of cross-section 
Flange width 
Depth of WF section 
Young’s modulus of elasticity 
Moment 
Moment applied at the end of a member 
= Moment at which yield point is reached in flexure 
Moment at which initial outer fiber yielding occurs when axial thrust 
is present 
Concentric axial thrust on a column 
Axial load corresponding to yield stress level over entire section 
Ratio of the maximum compressive residual stress, 7 rc, to the yield 
stress level, 7 
= Ratio of the maximum tension residual stress Ort to the yield stress 
level, 7 
= Section modulus 
Flange thickness 
Web thickness 
Depth of yield penetration in compression 
Depth of yield penetration in tension 
Distance from yield penetration in compression to extreme fiber in 
tension 
Plastic modulus 
= Ratio of the depth of yield penetration in compression to the total depth 


of the section (i.e., a= Yee) 
Ratio of the depth of yield penetration in tension to the total depth of 
the section (i.e., B 


= Centerline deflection of a beam-column loaded to produce single 
curvature 
Strain 
Strain corresponding to yield point stress 
Maximum residual compressive strain 
Maximum residual tensile strain 
Normal stress 
Lower yield point stress 
Maximum residual compressive stress 
Maximum residual tensile stress 
Variable normal stress on tension side of specimen (positive when in 
tension) 
Curvature at a section 
Curvature corresponding to initial outer fiber yielding (P = 0) 
Column erd rotation (at end A) 
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+ Table 4 
SUMMARY OF TEST CONDITIONS 
Test Loading P 
| Number Conditions Section ", r Type of Failure 
b 8 WF 31 0.49 55 Pending and Twist 
T-4 b 8 WF 31 0.12 55 Local Buckling 
ioe T-7 b 4WF 13 0.27 111 Pending and Twist 
T-9 b 4 WF 13 0.10 111 Bending and Twist 
T-10 b 4 WF 13 0.51 111 tending and Twist 
T-12 c 8 WF 31 0.12 55 ending and Twist 
iF T-13 d 8 WF 31 0.12 55 Local Buckling 
to T-14 a 8 WF 31 0.12 55 (Was not carried to 
failure) 
be T-16 c 8 WF 31 0.12 41 Eending, Twist and 
if 7 Local Fuckling 
1-17 b 4WF 13 0.12 55 Local Buckling 
oa T-19 © 8 WF 31 0.12 27 Pending and Twist 
' T-20 c 4WF 13 0.12 55 Bending and Twist 
T-21 b 4 WF 13 0.50 55 Local Puckling 
T-23 d 4 WF 13 0.12 84 Fending and Twist 
“ T-24 b 4WF 13 0.12 84 Pending, Twist and 
Local Buckling 
T-26 c 4 WF 13 0.12 84 Pending and Twist 
M,, 
M 
o 
P 
} * Slenderness Ratios Represent Free Lengths Between Knife-Edge 
or (Based on measured cross-section dimensions) 
330-20 
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Table 5 


Average Section Properties 


8 WE 31 


Measured 8.069 8.032 0.299 0.427 9.177 111.356 27.60 3.483 
Handbook 8.00 8.000 0.288 0.433 9.12 109.7 27.4 3.47 
Variation, % 0.86 0.40 3.82 1.39 0.63 1.51 0.73 0.39 


4 WF 13 


Measured 4.139 4.140 0.258 0.337 3.759 11.215 5.419 1.727 
Handbook 4.16 4.06 0.28 0.345 3.82 11.3 5.45 1.72 
Variation, % 5.05 4.93 7.86 2.32 1.59 0.75 0.57 0.41 


Material Properties 


Section 


8 WF 31 40.0 ksi 29,900 ksi 
4W 13 39.5 ksi 200 ksi 
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Strain Diagram Stress Diagram 
(a) (b) 


ig.! Elastic Stress and Strain Distribution 


Fig.2 Idealized Stress-Strain Diagram 


Strain Diagrom Stress Diagram 
(a) (b) 
Fig. 3 Stress and Strain Distribution Patterns 
when Strain Exceeds the Elastic Limit 
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Comp. | Tension 


Fig. 4 Hypothetical Cooling Residual Stress Pattern 
for a Flat Piate 


Fig. 5 Typical Cooling Residual Stress Pattern 
for a W Section 
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Fig.6 Measured Residual Stress Pattern 
for the 8W Section 


Fig. 7 Assumed Cooling Residual Stress Patterns 
for the Wide Flange Section 
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Original 


Cold Bend 


Final 


Fig.8 Diagrammatic Representation of the 
"Gag" used in Straightening Rolled W Shapes 


Tension | Compression 


Fig.9 Residual Stress Pattern due to 
Cold Bending 
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Residval Strain Applied Strain Resulting Strain 
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Residual Stress Kdpited Stress Resulting Stress 
(a) (b) (c) 


Fig.10 Influence of Cooling Residual Stresses on 
Elastic Stress and Strain Distribution 
(Plate Sections ) 


Residual Stress Applied Stress Resulting Stress 
(a) (b) (c) 


Fig.!| Influence of Cooling Residual Stresses on 
Stress and Strain Distribution when Strains 
Exceed the Elastic Limit (Plate Sections) 
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Tension Compression T,C 


Top Flange 
(Compression) 


Bottom Flange 
(Tension) 


Comp. Flange Partial Yield Complete Yield | Complete Yield 
Tension Flange § Below Yield Pt. | At Yield Point | Partial Yield 
4. 


l. 2. 3. 


Note: Applied Stress Patterns ore Cross Hatched 


Fig.12 Influence of Cooling Residual Stress on Stress 
Distribution for W Sections (Strong 


Axis Bending) 
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Top Flange 
(Compression) 


Bottom Flange 
(Tension) 


Comp. Flange At Yield Point Partially Yielded | Complete Yield 
Tension Flange} Below Yield Point| At Yield Point Partially Yielded 


Note: Applied Stress Potterns are Cross Hatched 


Fig. 13 Influence of Cold Bending Residual Stress 
on Stress Distribution for W Sections 


(Strong Axis Bending) 
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Strain Diagram Stress Diagram 


Fig. 14 Elastic Stress ond Strain Diagrams 


Pee 


Strain Diagram Stress Diagram 


Fig. 15 Stress and Strain Diagrams 
when Strains Exceed the Elastic 
Limit 
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Fig. 16 Typical Progressive Yield Conditions 
Due to Beam- Column Loading 
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Fig.17 Sample Set of Assumed Yield Penetration 
Conditions for Determination of Auxiliary 
M-¢ and P-¢ Curves 
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Fig. 18 Auxiliary M-¢ and P-¢ Curves for an Assumed 
Yield Penetration of 1/4 the Depth of the Section 
(8W 31) 
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Fig.19 Typicol Set of M-¢ and P-¢"Auxiliory Curves" 
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Fig. 21 Enlarged M-¢ and P-¢ “Auxiliary Curves" 
(8W 31) 
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Fig 22 M-¢ Curves for the 8W 3! Section 
(Strong Axis Bending) 
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Fig.23 “Auxiliory Curves” Including the infiuence of 
Residual Stress 
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* See Ref.9 For Defination of Strength Terms 
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